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Wird  die  Rotation  eines  festen  Körpers  durch  einen 
Körper  gestört,  der  sich  in  einer  Ebene  bewegt,  so  wird  diese 
Ebene  zweckmässig  aJs  XF-Ebene  gewählt.  Der  Ausdruck 
fiir  das  Potential  vereinfacht  sich  dabei,  indem  alle  mit  dem 
Sinus  der  Breite  des  störenden  Körpers  multiplizierten  Glieder 
verschwinden.  Diese  Vereinfachung  ist  erlaubt  so  of  t  es  sich 
um  die  Störungen  der  Sonne  auf  die  Rotation  der  Planeten 
handelt,  und  ich  werde  deswegen  dies  Problem  kurz  als  das 
planetarische  Roiationsjrrohlem  bezeichnen. 

Das  planetarische  Rotationsproblem  besitzt  einige  parti- 
kuliäre  Integrale,  die  ich  hier  untersuchen  will.  Ich  bediene 
mich  der  Veränderlichen  ij,  ^^,  u,  v,  p,  q,  die  ich  an  anderer 
Stelle  definiere.  ^  Die  Störungsfunktion  ist  nach  Potenzen 
von  u,  v,  p,  q  entwickelt.  Wird  mit  P,  {u,  v,  p,  q)  öder  kurz 
Pr  eine  Potenzreihe  in  u,  v,  p,  q  bezeichnet,  in  welcher  die 
Summe  der  Exponenten  der  niedrigsten  Dimensionen  in  der 
Potenzreihe   gleich  )•  ist,    so  känn   man   die   charakteristische 

'  »Entwickelung  des  Potentials  im  Rotationsproblem.i  Ich  werde  iin 
Folgenden  diesen  Aufsatz  mit  Meddel.  n:o  33  bezeichnen. 

Aikiv  för  matematik,  astronomi  och  fysik.     Bd  4.     N:o  12.  1 
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Funktion  H'  in  der  Form 


3  u 


Ä"'  =  20  "  2  ?  1^^^  "^  ^^^^  '^'  +  hB  sin»  ?y  J  — 
—  ii^  +  P^  {u,  v,  p,  q) 

schreiben. 

Die  Differentialgleichungen  der  Bewegung  sind    (Medc 
N:o  33) 

dt  ~  r7ijj'     (f^  ~       d^^' 

du_<2H^      dv^__<lH[ 
dt        dv'     dt  ~       du ' 

dp^OH^      ^  _  _  'W 

dt        äq  '      dt  <)p ' 


(2) 


Es  erhellt    aus   diesen  Gleichungen   und   aus   der  obi^ 
Form   fiir  H',    dass   das  planetarische   Rotationsproblem 
partikularen  Lösungen 


(3) 


M  =  v  =  p^<7  =  0 


hat,  wenn  nämlich  ij  und  r^^  aus  den  ersten  Gleichungen 
bestimmt  werden.     Wir  setzen  zur  Abkiirzung 


(4)         i^o  = 
so  dass 


W~2%'^^^"^  ^^®^ ''»  -^^2  B  sin^ ».,]  —  ^-,  j^ , 


ist. 


(5) 


Dann  ist 


H'  =  Fo  +  P,  (ti,  v,  p,  q) 


Fiihrt    man    die    Differentiationen    aus,    so    lauten    di 
Gleichungen,  wenn  man  beachtet,  dass 


k,A—k^B  =  B  —  A 


ist, 
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(5*) 


2  r^ 


di.       3  u  ,-,        ..  ^ 

7^  =  9  *3  (^  —  ^)  2  cos  »;,  sm  r^^ 

dt 


i,      dl 
C^df 


Die  Länge  A  und  der  Radius  Vector  r  des  störenden 
Planeten  sind  gewisse  Funktionen  der  Zeit,  die  wir  im  all- 
gemeinen  als  bekannt  voraussetzen  können.  Wir  fangen  mit 
dem  einfachsten  Fall  an,  dass  nämlich  der  störende  Körper 
sich  in  einem  Kreis  um  den  rotierenden  Körper  bewegt. 
Bezeichnen  wir  die  mittlere  Bewegung  des  störenden  Körpers 
mit  n,  so  hat  man  dann 


(6) 


dl 
dt 


u 


u  +  «i 


wo   III   di®  Masse   des   rotierenden  Körpers  bezeichnet.     Wir 
setzen 


(7) 


3  ^ 


u  B 


C 


Die  Grösse  r  ist  also  eine  Konstante,  wenn  der  störende 
Körper  sich  in  einem  Kreis  bewegt.  Die  Gleichungen  (5*) 
lauten  dann 


(8) 


~  =  6  V  cos  »^,  sm  r]i , 


dr^i 
[  dt 


+  n. 


Hier    findet   man    unmittelbar    die    beiden    partikularen 
Lösungen 


\ 


und     /,,  =  ±  90°,  il  =  Cn. 

Beivegt  sich   der  störende  Körper  in  einem  Kreis,    hat  also 
■s   pla7ietafische    Rotations  problem   die   partikularen    Integrale 
v  =  p  =  q  =  O,  ^1  =  Cn  und  r^i  ==  O  öder  »ji  =  ±  90°. 
Die  mechanische  Bedeutung   dieser  Lösungen  sind  leicht 
zu  finden.     Da   w  =  v  =  p  =  g  =  O  ist,    so   ist  0^  und  e  auch 
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gleich  Null,  so  dass  die  Rotationsachse  mit  der  Achse  des  Tråg- 
heitsmomentes  C  zusammen fällt. 

Was  i]i  betrifft,  so  hat  man  nach  der  Definition  (Meddel.  33) 

1]  =  U^  +  U2  —  U^  +  K . 

Wenn  die  drei  Achsen  Z,  z  und  Z^  zusammenfallen,  so 
ist  (Meddel.  31,  Fig.  1) 

^,==180°,     w,  — ^t3=  180°  + a;Z, 

so  dass 

i^i  =  A  —  Xx 

Die  Grösse  *;i  bezeichnet  also  hier  den  Unterschied 
zwischen  der  Länge  der  Planeten  und  der  Länge  der  Träg- 
heitsachse  x  in  der  Z7-Ebene.  (Die  Trägheitsachse  x  ent- 
spricht  dem  Trägheitsmoment  A). 

Die  beiden  Lösungen  »^,  =  O  und  /;i  =  90°  bedeuten  also 
mechanisch,  dass  die  kleinste  bezw.  die  mittlere  Trägheits- 
achse gegen  den  störenden  Körper  gerichtet  sind.  Man  hat 
es  mit  sog.  gebundener  Rotation  zu  tun. 

Diese  partikularen  Integrale  sind  indessen  nur  Spezial- 
fälle  der  allgemeinen  Lösungen  der  Gleichungen  (8),  die  leicht 
zu  finden  sind.  Wird  nämlich  die  zweite  Gleichung  (8) 
differentiiert,    so  erhält  man   fiir  *yi  die  Differentialgleichung 


Diese  gibt 


d^  >/,  , 

-Jt2  =—  '"  co«  %  sin  »il 


dt 


wo  E  eine  Integrationskonstante  bezeichnet. 

Das  Integral  dieser  Gleichung  ist  verschieden  je  nachdem 
der  Wert  von  E  grösser  öder  kleiner  als  >'-  ist.     Ist 

1)  E>  r^ 

so  ist  die  Lösung 

(9)  rj,  ^ami^^Et^  y)-,  Mod.  J    . 

'  E 
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Ist  dagegen 


i 


2)         E  <  i>\ 
erhält  man  eine  reelle  Lösung  von  der  Form 


^  E  ^  F 

( 10)  sin  »;,  =  —  sn{vt-\-  y);  Mod.  — -  • 

Im  letzteren  Falle  haben  wir  also  Libration  in  »ji  um  den 
Wert  Null.     Der  Maximalwert  von  /^j   ist  durch  die  Relation 


sm  r 


Ve 


/i 


gegeben.  Die  Achse  des  kleinsten  Trägheitsmomentes  schwankt 
periodisch  um  den  Radius  Vector  des  störenden  Körpers 
ganz  in  derselben  Weise  —  und  nach  ähnlichen  Gesetzen  — 
wie  das  einfache  Pendel  um  die  Lotlinie.  Beispiele  in  der 
Natur  dieser  Art  sind  die  Libration  des  Mondes  und,  nach 
ScHiAPARELLi,  die  Rotation  des  Mercurs  und  der  Venus. 

Ist  E  =  v^  so  findet  Limitation  statt.  Man  hat  in  der 
Tat  dann,  wie  aus  (9)  öder  (10)  hervorgeht, 

(^eht  t  gegen  +  c>o,  geht  sin  »;i  gegen  die  Einheit,  und  also 
*;i  gegen  90°.  Die  mittlere  Trägheitsachse  nähert  sich  asympto- 
tisch  dem  Radius  Vector  des  störenden  Körpers.  Dies  ent- 
spricht  der  partikularen  Lösung  »;i  =  90°. 

Es  erhellt  aus  den  obigen  Ausdriicken  unmittelbar,  dass 
die  Libration  um  ?;i  ==  O  stabil,  um  jj^  =  90°  dagegen  unstabil 
ist. 

Ist   E  >  v^   so   wächst    »^j    mit    der  Zeit    ins  Unendliche. 

Wird  //i  in  eine  FouRiER'sche  Reihe  entwickelt,  so  erhält 
man  fiir 


1)  E>  v' 


.r(l^j^«  +  7)  ,   ^V'l      q'  ^^^^Jt}fEt±y)^ 


2R  ^"^5l  +  g2.-"^  It 


6      AEKIV  FÖR  MATEMATIK,  ASTRONOMI  O.  FYSIK.      BD  4.    N:0  12. 

Ist  dagegen 

2)  E<  v^ 


so  erhält  man 


und  hinaus 


^=^yEcn{vt  +  y) 


_4y   _! VJ^-i     ,    {2s  —  \)7t(vt  +  y) 

Ein  charakteristischer  Unterschied  zwischen  den  beiden 
Lösungen  verdient  besonders  hervorgehoben  zu  werden.  Er 
betrifft  die  Grösse  der  Periode  der  Schwankungen.  Ist 
nämlich 

1)  E>  v^ 

so  erhält  die  Schwankungsperiode  den  Wert 

2K 

Ve 

und  fiir 

2)  E  <  v^ 
den  Wert 

4Z 
v 

Fiir  sehr  hohe  Werte  der  Integrationskonstante  E  hat  mit- 
hin  die  Schwankungsperiode  genähert  den  Wert 

yr 

Vi 

urid  fiir  sehr  kleine  Werte  dieser  Integratiönskonstante  hat  sie 
genähert  den  Wert 

v 

Im  vorigen  Ealle  nähert  sich  die  Schwankungsperiode 
dem  Werte  NuU,  im  letzteren  Falle  dem  Werte  Inw,  der 
vom|]  Unterschied   der  Trägheitsmomente  B  und  A  abhängig 
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ist.     Diesen   Grenzwert  der  Schwankungsperiode   können  wir 
in  der  Form 


k 


ary^i^ 


schreiben,    wenn  2T  die  Umlaufszeit   des   störenden  Körpers 
bezeichnet. 

Der  Minimalwert  von  {B  —  A):  C  ist  Niill  (Rotations- 
körper  um  die  grösste  Trägheitsachse  rotierend),  der  Maximal- 
wert  ist  Eins  (fadenförmiger  Körper).  Die  Schwankungsperiode 
ist  also  grösser  als 

^  =  0.577X2T. 

Nachdem  ?;i  gefunden  worden  ist,  erhält  man  ^j  aus  einer 
der  Gleichungen  (8).     Aus 

dB 

-~  =  C  r^  cos  >ji  sin  1^1 

erhält  man  fiir 

1)  E>i'^ 

^  =  C  r^cni^E t  +  y)  sn  {^Et  +  y) 
dt 

also 

i-j  =  —  O  ^Edn{^Et  +  y)  +  Konstante, 

wo  die  Konstante  nach  der  zweiten  Gleichung  (8)  den  Wert 
■w  hat.     Ist 

2)  E<v^ 
\o  erhält  man 

^  =  Cv  ^Edn  {vt  +  y)  sn  {vt  +  y) 
dt 

und 

§'j  =  —  C  ^E  en  {v  ti-  y)  +  Konstante, 

70  die  Konstante  wieder  den  Wert  Cn  hat. 
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Wir  haben  bis  jetzt  angenommeii,  dass  der  störende 
Körper  sich  um  den  rotierenden  Körper  in  einem  Kreis  bewegt. 
Ist  dies  nicht  der  Fall,  so  kompliziert  sich  die  Aufgabe.  Die 
Differentialgleicliung  fiir  r.^  hat  dann  die  Form 

(11)  ^  +  >'2  C08  ij,  sm  j^i  =  ^  , 

lind  hier  sind  sowohl  r-  wie  I  Funktionen  der  Zeit. 

Indem  ich  in  der  ersten  Linie  den  Librationsfall  vor 
Augen  habe,  so  ist  ersichtHch,  dass  die  Veränderungen  von 
v\  welche  von  Veränderungen  im  Radius  Vector  herriihren, 
von  geringerem  Einfluss  sind,  als  die  Veränderungen  in  A. 
Die  Veränderungen  in  h'^  sind  nämlich  mit  der  als  klein 
angenommenen  Grösse  sin  rji  multipliziert.  Wird  fiir  A  die 
Form 

I  =  nt  +  Co  + 2  ÄGOS  {ht  +  c) 

ängen ommen,  so  känn  man  in  vielen  Fallen  die  Differential- 
gleichung  (11)  auf  die  lineare  Form 

bringen,  und  das  Integral  lautet  dann 

(12)  iq^  ==  Ao  COS{vt  +  So)+2iJ)2__  ^2  G08{bt  +  C), 

vvo  A,y  und  Bq  die  Integrationskonstanten  bezeichnen. 

Diese  Formel  gibt  eine  gute  Annäherung  zum  gesuchten 
Integral  so  oft  keine  von  den  Grössen  b  genähert  gleich  y 
sind.  In  letzterem  Fall  entstehen  aber  Schwierigkeiten,  die 
man  bis  jetzt  nicht  geniigend  untersucht  hat.  Die  Unter- 
suchungen,  die  mir  bekannt  sind,  geben  wenigstens  keinen 
Aufschluss  iiber  das  Verhalten  des  Integrales  fiir  solche 
Werte  des  Parameters  r^,  die  in  der  Nähe  des  Wertes  p^  =  b^ 
liegen. 

Ich  nehme  an,  dass  es  sich  um  die  Differentialgleichung 

(j     fl 

(13)  -yy^-  +  v-  cos  *^  1  sin  *^,  --=  —  b'^A  cos  b  t 
handelt,  so  dass  nur  ein  Glied  in  'A  mitgenommen  wird. 
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Zu  dieser  Gleichung  existiert  immer  ein  partikulares 
Integral  mit  der  Periode  2 .c :  b.  Ich  will  dies  Integral  näher 
in  Betracht  ziehen.     Ich  setze  also 

( 13*)  /;i  =  «i  cos  &^  +  «3  cos  3  &<  +  «5  cos  5bt  -] 

und   suche    die    Koeffizienten   a^,  a^,  a.^   u.  s.  w.    durch    Re- 
kursionsformeln  zu  bestimmen.     Man  hat 


und 


1    .  2  2 

cos  /;  sm  r^,  =  g  »m  2  >^,  =  i^,  —  3^1  +  j^  vl 


t^l=  .aliScosbt  +  cosSbt) 

3 

+  ^«!f'3  {cos  bt  +  2  cos  3  6^  +  cos5&0 

3 

+  ,  «i«3  (2cos&^  +  cos  5bt  +  cos  7 60 

+     «^(3cos3  6«  +  cos9&<) 


Hieraus  entstehen  die  folgenden  Rekursionsformeln : 

—  b' A  -=  —  6-«i  +  r2(«i  —  2 «i  —  2  "1  "3  —  «i «3  H ) 


r 


0  =  — 9&2«3+  i,^c(s  —  ^ccl  —  alc(s  —  ~^al  +  ■■■) 


Es  ist  wohl  kaum  zu  erwarten,  dass  man  eine  strenge 
Auflösung  dieser  Gleichungen  finden  könnte.  Um  eine  Vor- 
stellung  von  der  möglichen  Form  des  Integrales  zu  erhalten^ 
bemerke  ich,  dass  nach  der  zweiten  Gleichung  a^,  fiir  i^^  =  ö^, 
den  genäherten  Wert 


nat,  ^ 


«i_ 


at7  woraus  man  findet,  dass  «3  sehr  klein  ist,  wenn  «,  einen 
kleinen  Wert  hat.    Setzt  man  dementsprechend  in  der  ersten 
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Gleichung  a^  =  O,  so  erhält  man  fiir  «i  die  Gleichung 
(14)  (-,-l)«,  +  2««  =  -,^. 

Wäre  v^=b-,  so  wiirde  man  hieraus 

erhalten . 

Dies  ist  aber  nicht  der  grösste  Wert  von  a^,  der  auf- 
treten  känn.    Diesen  trifft  man  vielmehr,  wenn  die  Gleichung 

(14)  eine  Doppelwurzel  besitzt.     Schreibt  man 

so  ist 

(15)  al  —  2{l  —  z)ai  —  2zA==0 

und  diese  Gleichung  hat  eine  Doppelwurzel 

8  

«!  =  —  ^zA, 
wenn  z  die  Gleichung 

(16)  z'Ä^-  =  §j{l-^)' 

befriedigt. 

Lässt  man  z  von  Null  bis  oo  variieren,  so  findet  man, 
dass  diese  Gleichung  eine  einzige  positive  Wurzel  besitzt. 
Diese  Wurzel  ist  zwischen  z  =  O  und  z  =  +  1  gelegen.  Setzt 
man  nämlich 

f{z)==§j{z-ir  +  z'A', 

so  findet  man,  dass  /(O)  negativ,  /(I)  positiv  und  f  (z)  zwischen 
2  =  0  und  2  =  +  c\3  positiv  ist.  Es  gibt  also  eine  einzige 
positive  Wurzel.  Diese  Wurzel  —  die  wir  mit  Zq  bezeichnen  — 
lässt  sich  nach  Potenzen  von  A^^i^  entwickeln  und  zwar 
bekommt  man 

z^  =  l-lA'h  +  A'l^---- 


I 
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Ist  A  klein,  so  ist  Zq  nahe  gleich  der  Einheit.  Der  eiit- 
sprechende  Wert  von  «i  ist 

welche  Wurzel  doppelt  ist.     Ausserdem  besitzt  (15)   die  ein- 
fache  Wurzel 

3   

ai  =  +  2^Zo^- 

Ist  z>Zq  hat  (15)  eine  einzige,  und  zwar  positive,  Wurzel, 
die  sich  fiir  wachsende  z  dem  Werte  «i  =  A  nähert.  Ist 
z<Zo  so  hat  (15)  3  reelle  Wurzeln,  von  denen  eine  positiv, 
die  beiden  anderen  negativ  sind.  Fiir  2  =  0  hat  man  die 
Wurzeln  «!  =  O  und  a^  =  ±  ^2 . 

Die  periodische  Lösung  (13*) 
hat  •  also  die  merkwiirdige 
Eigenschaft  eine  diskontinuier- 
liche  Funktion  des  Parameters 
z  zu  sein.  Fangen  wir  näm- 
lich  mit  dem  Werte  s  =  O  an, 
so  hat  (13)  die  Lösung  r/i  =  0. 
Wenn  z  wächst,  so  wachsen 
die  Koeffizienten  (wenigstens 
gilt  dies  von  dem  ersten  Glied 
in  der  Reihe  (13*))  in  der 
periodischen  Lösung,  bis  wir 
fiir  z  =  Zo  die  Lösung 

s 

>^j  =  —^Zf,A  cos  6<  +  «3  cos  3  &<  +  •  •  • 

erhalten.     Wenn   nun  z>Zo  wird,    so   haben  wir  es    mit   der 
periodischen  Lösung 

3  

r^i  =  2^ZoA  cos  bt  +  asCOsSbt  +  •  •  • 

zu  tun,    d.  h.  die  Amplitude    wird   verdoppelt   und  wechselt 
das  Zeichen. 

Die  hier  behandelten  Lösungen  haben  ihre  Anwendung 
bei  dem  Studium  der  Bewegung  eines  einfachen  Pendels,  das 
Störungen  unterworfen  ist,  welche  nahe  dieselbe  Periode 
besitzen  wie  die  freien  Schwingungen  des  Pendels.     Die  hier 
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angewandte  Methode  ist  zwar  in  mathematischer  Hinsicht 
sehr  mangelhaft,  diirfte  aber  eine  genäherte  Aufgabe  der 
Lösung  geben. 

Es  ist  fiir  die  Theorie  der  Libration  des  Mondes  von 
Interesse,  die  numerischen  Beträge  der  hier  behandelten 
Störungen  zu  untersuchen.  Das  Glied  in  der  Mondlänge,  das 
hier  in  Betracht  kommt,  ist  die  sog.  jährliche  Gleichung, 
welche  nach  Delaunay  die  Form 

—  668'.' 9  sin  O 

hat.  Man  hat  also  A  =  —  668'.'9 ,  öder  im  Radius  ausgedriickt 
A= — 0.00324.  Wenn  die  Trägheitsmomente  des  Mondes 
solche  Werte  liaben  wiirden,  dass  die  Periode  der  physischen 
Libration  genau  ein  Jahr  betragen  wiirde,  so  wiirde  man  in 
r]i  ein  Glied  erhalten  vom  Betrage 

—  11! 9 9  sin  O. 

Im  singulären  Punkt  Zg  springt  der  Koeffizient  dieses 
Gliedes  von   +  8!4  8  zu  —  16^96  iiber. 

Da  die  Beobachtungen,  die  man  zu  verschiedenen  Zeiten 
angestellt  hat,  um  den  Wert  des  GHedes  im  Ausdrucke  fiir 
die  Libration  des  Mondes,  das  die  Periode  eines  Jahres  hat, 
zu  bestimmen,  insofern  iibereinstimmen,  dass  der  Koeffizient 
dieses  Gliedes  höchstens  einige  Bogenminuten  betragen  känn, 
so  scheint  es  also  ausgeschlossen  zu  sein,  dass  die  Periode  der 
physischen  Libration  die  Länge  von  einem  Jahr  hahen  könnte. 

Wir  haben  bis  jetzt  nur  die  partikulare  periodische 
Lösung  der  Gleichung  (13)  betrachtet.  Um  die  allgemeine 
Lösung  zu  erhalten,  miissen  wir  noch  die  entsprechenden 
Variationsgleichungen  in  Betracht  ziehen.     Setzt  man 

tj  =  tj"^  +  Jr] 

und  verstehen  unter  j^"  die  partikulare  Lösung  von  (13),  die 
eine  Periode  gleich  2. t -.b  hat,  so  erhält  die  Gleichung  fiir 
Jtj  die  Form 

^+^^cos2,,;^,;  =  0. 

Diese  gehört  zu  den  linearen  Differentialgleichungen  mit 
periodischen  Koeffizienteii,  und  das  Integral  ist  von  der  Form 
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(17)  Jr^^  ^Ugcoslw  +  2sbt], 


i 


3VO 

w  =  r  {]  —  o )  /  +  .T 


ist  lind  die  Grösse  a  genähert  den  Wert 

hat  (man  vergleiche  meine  Vorlesungen  iiber  die  »Mechanik 
des  Himmels»  S.  41).  Wenn  r-  nahe  gleich  b-  ist,  so  ist 
dieser  Ausdruck  nicht  mebr  giiltig. 

Fiigt  man  zu  (17)  das  partikulare  Integral  (13*)  hinzu, 
so  hat  man  die  allgemeine  Lösung  der  Differentialgleichung  (13). 

Die  partikularen  Integrale,  die  wir  hier  betrachtet  haben, 
werden  zweckmässig  als  Ausgangslösungen  benutzt,  wenn  man 
die  Rotation  des  Mondes  untersuchen  will.  Ich  habe  bei  der 
Untersuchung  dieses  Problems  diesen  Weg  eingeschlagen  und 
dadurch  eine  grosse  Einfachheit  und  tlbersichtlichkeit  in  der 
Behandlung  des  Problems  erreicht. 

Ist  der  rotierende  Körper  ein  Rotationskörper,  so  dass 
A  =-B  ist,  so  besitzt  das  planetarische  Rotationsproblem  — 
vorausgesetzt,  dass  der  störende  Körper  sich  in  einem  Kreise 
bewegt  —  ausser  dem  Integrale  H'  ==  Constans  noch  ein  Inte- 
gral (nämlich  »i  =  Constans).  Die  Differentialgleichungen 
können  dann  auf  zwei  Freiheitsgrade  reduziert  werden.  Dies 
findet  seine  Anwendung  in  der  Theorie  der  Rotation  der  Erde. 
Das  Integral  cq  =  Constans  erlaubt  wichtige  Schliisse  iiber 
die  Stabilität  der  Rotation  der  Erde  zu  ziehen. 


Tryckt  den  28  januari  1908. 

Uppsala  190S.    Almqvist  &  WikselU  Boktryckeri-.\.-B. 
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